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Abstract 

The physics underlying the magnetization process of quantum antiferromagnets is revisited from the 
viewpoint of geometric phases. A continuum variant of the Lieb-Schultz-Mattis-type approach to the prob- 
lem is put forth, where the commensurability condition of Oshikawa et al derives from a Berry connection 
formulation of the system's crystal momentum. We then go on to formulate an effective field theory which 
can deal with higher dimensional cases as well. We find that a topological term, whose principle function is 
to assign Berry phase factors to space-time vortex objects, ultimately controls the magnetic behavior of the 
system. We further show how our effective action maps into a Z2 gauge theory under certain conditions, 
which in turn allows for the occurrence of a fractionalized phase with topological order. 

PACS numbers: 75.10.Jm 
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I. INTRODUCTION 



Antiferromagnets subjected to an external magnetic field have attracted considerable attention 
over the years, largely owing to interesting features which they share with a rather broad class of 
quantum many-body systems situated on a lattice. The series of quantum phase transitions en- 
countered as the field strength is varied may be viewed as a spin-analogue of the superfluid-Mott 

nnr 

insulator transitions observed in boson Hubbard models [11121130, which model e.g. Josephson junc- 
tion arrays[|4|] and optical lattices[5]. Magnetization plateaus emerging at simple fractions of the 
saturated magnetization are reminiscent of the quantum Hall effects |6i LZD - Recent developments 
reveal the subject to be interconnected with an even richer variety of issues such as the Luttinger 
theorem[8], electric polarization in solids [uj], and fractionalization of quantum numbers in 



dimensions greater than onellll 



In this article we revisit this problem from yet another perspective, i.e. that of geometric 
phases II 13l 1 1411. The basis of our arguments rests only on rather general properties of Berry con- 
nections, topological terms and boson-vortex duality, much of which are valid in any dimension. 
We aim to shed new light on the above body of interrelated phenomena, with particular empha- 
sis on providing a workable format with which to pursue further exotica in quantum spin liquids, 
anticipated to emerge along this line of study. 



H. OUTLINE OF THIS PAPER 



In view of the fairly technical nature of the presentation to follow, we wish to highlight in this 
section its main contents, placing them in context with previous work. 

Section Hn] is devoted to a geometrical reinterpretation of the well-known d=l result of Os- 
hikawa, Yamanaka and Affleck (OYA)|f|], which yields a quantum- mechanical constraint on the 
possible values which the magnetization can assume at plateaus, i.e. in a spin-gapped regime. 
(In this article the notation d is reserved to denote the spatial dimensionality.) The OYA work 
(reviewed in subsection IIII Al) builds on arguments initiated by Lieb, Schulz and Mattisl 15. [l^l 
(LSM), wherein the central step consists of comparing the crystal momenta of the ground state and 
a candidate low-lying excited state, constructed from the ground state by applying a slow twist to 
the spins. (The resulting quantization rule agrees with that obtained by bosonization methods 
QJ/ZD -) Meanwhile there is a simple geometric formulation of the crystal momentum of aferro- 
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magnetic spin chain due to Haldane[18], which incorporates the language of spin Berry phases. 
We show that this framework can be adapted to our problem involving antiferromagnets, provided 
the spin moments are partially polarized due to the magnetic field. Carrying out the LSM proce- 
dure in this geometric language then simply amounts to evaluating the difference between Berry 
phases associated with untwisted and twisted spin configurations. We find that this indeed leads 
immediately to the celebrated OYA rule. 

In section |IV] we show that a similar geometric structure underlies the magnetization properties 
of antiferromagnets in arbitrary spatial dimensions. Here there are no obvious substitutes for the 
LSM scheme or the bosonization technique (both of which are by design specific to d=l), and the 



generalization of the OYA result to spatial dimensions larger than unity II 1911 is known to be a subtle 



problem. Our strategy is to build instead on the semiclassical picture of the previous subsection 
and to derive a low-energy effective theory for partially polarized antiferromagnets. With the 
magnitude of the spin polarization essentially fixed by the magnetic field, the long wavelength 
physics mainly involves the orientational fluctuation of the residual planar staggered moment, 
lying in the plane perpendicular to the field. Symmetry thus dictates that our effective action 
should be a variant of the quantum XY model. A careful derivation confirms this expectation; we 
find however that the imaginary-time action for this XY model also contains a purely imaginary 
topological term (proportional to the time-derivative of the phase variable) whose coefficient is 
given by S — m, with S the spin quantum number and m the magnetization per site. This perhaps 
is not surprising if one recalls that a complex- valued low energy action also arises out of the 



20, 



210. 



closely-related boson-Hubbard model away from commensurate filling factorsyj, 

The later half of section [IV] thus explores the consequence of complexifying our action by the 
addition of such a term to the effective XY model. We know from earlier work on the hydrodynam- 
ical properties of superfluids and superconductors (where actions of the same form appear) that 
the topological term crucially influences the quantum dynamics of phase vortices by subjecting 



them to a Magnus -type force |22L|23 



24J]. At the heart of this phenomenon is an Aharonov-Bohm 



(AB)-like quantum phase interference which renders each vortex to act like a charged particle im- 
mersed in a fictitious magnetic field whose strength is fixed by the topological term. (The Magnus 
force is then simply understood to be a pseudo-Lorentz force.) Since vortex events (point vortices 
in (l+l)d, vortex loops in (2+l)d, and vortex sheets in (3+l)d) are the most relevant disordering 
agents of an XY model, it is clearly important to see how such an interference effect manifests 
itself in the present problem. For this purpose we incorporate simple duality techniques which 



3 



enables us to visualize the physics directly in the language of vortices. (As later explained, the 
approaches taken here as well as in section |V] are intimately related to the recent work of Balents 
et al[25] who (in a context slightly different from ours) run a detailed projective symmetry group 
analysis to classify the possible ground states which the vortices of the XY model can assume.) 
We find that the pseudo-AB effect generates Berry phase factors associated with the vortex events, 
which in turn governs (given a specified valued of S — m) whether the condensation of these 
topological defects is allowed as a result of constructive interference, or alternatively, a destruc- 
tive interference renders the vortex configurations irrelevant. (The former case will imply spin 
disordering and hence the formation of a magnetization plateau.) It is easy to recognize that the 
former can happen for S — m G Z, the condition under which all vortex Berry phases become 
trivial. Defects of arbitrary vorticity (most importantly singly quantized vortices) will then be able 
to condense. A more subtle situation arises when the value of S — m is a rational number, in which 
case vortices must be multiply quantized in order to condense. With this physical picture (which 
links Berry phases to the magnetization behavior) in hand, it is also interesting to study the effect 
of Dzyaloshinskii-Moriya type perturbations which can modify the Berry phase factors. A brief 
discussion on this issue is given at the end of this section. 

In section V we look into a particularly important class of problems where the quantity S — m 
is a half-odd integer, i.e. S — m G Z + 1/2. The central observation to make here is that the 
system can dynamically acquire a Z 2 gauge symmetry. We thus come into contact with the work 
of Sachdev and Park[26] who, in the course of their search for novel spin liquids in the absence 
of a magnetic field arrive at the same effective theory (for easy-plane antiferromagnets) - a lattice 
XY model coupled to an Ising gauge theory with a Z 2 -valued Berry phase term. This theory has a 



2711 . and we discuss how this arises 



phase diagram which accommodates a fractionalized phase[26, 
within our magnetization problem. Interestingly, we find that this derivation generalizes to the case 
where S — mis set at other simple rational numbers; for instance when S — m G Z + 1/3, we are 
lead to a similar effective theory coupled to a Z 3 gauge field with a Z 3 -valued Berry phase term. 
The results of this section thus illustrates the unexpectedly rich phase structure of square lattice 
antiferromagnets in a magnetic field, suggesting them to be a promising place to seek exotic spin 
liquid states. 

Appendix lAl supplements section |In] while the details of the duality methods used in section ITVl 
can be found in Appendix |Bl 
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III. GEOMETRICAL INTERPRETATION OF THE LIEB-SCHULTZ-MATTIS ARGUMENT 



A. A brief summary of Lieb-Schultz-Mattis 

The LSM argument! 15. gives us an important insight into how low-lying excitations in 
one dimensional systems are constructed. Let us consider a translationally invariant d=l quantum 
system with a unique ground state |G.S.). The key idea of LSM is that the following twisted state 

C/ LSM |G.S.> = e^^|G.S.) (1) 

is made up of low-lying excited states (provided that the crystal momentum satisfies a certain 
condition which will be fixed shortly), and should hence be a useful reference state for extracting 
information on the low energy spectrum of the system. To evaluate the momentum shift caused by 
t^LSM, it is convenient to consider the phase acquired in the process of executing a 1-site translation 

f (f S n f - 1 = S„_i). Since 

TUlsmT- 1 = e-^^-^sM , (2) 
one can readily see that ?7lsm shifts the crystal momentum P by 

5P^ M = -—Y,(S-S^) = -2ir{S-m) mod27r, (3) 
i=i 

where we have introduced the magnetization (density) m = V - Sj/N. The above expression 
implies that if S — m ^ Z the twisted state £/lsm|G.S.) is orthogonal to the ground state |G.S.), 
that is, the twisted state consists of excited states of the finite-size system. Moreover, it is not 
difficult to show that for generic Hamiltonians with short-range interactions the energy of the 
twisted state decreases like 1/N as the system size grows N /* oo. Thus we are lead to an 
interesting dichotomy; if S — m ^ Z, the system in the thermodynamic limit either has gapless 
excitations over the unique ground state or has several degenerate ground states corresponding 
to a spontaneous breaking of (translational) symmetry. (The latter possibility becomes relevant 
when S — m = p/q (where integers p and q are coprime), as one easily sees by repeating the 
above twist operation q times.) This, combined with a complimentary commensurability argument 
based on bosonization, yields the so-called quantization condition of magnetization plateaus^: 
Qg.s.(S — m) G Z (Qg.s.: period of the infinite-size ground state). 

In a strict sense, the above LSM argument alone tells us nothing about the presence/absence of 
plateaus as low-lying excitations created by the LSM twist do not change the total magnetization 



Sj. Moreover, although a similar conclusion has been reached[19] in higher dimensions as 



well, it is not totally clear why the same sort of quantization condition exists regardless of di- 
mensionality while the LSM argument is valid only in one dimension. In an attempt to clarify 
such issues, and to set the stage for unearthing further exotic properties, we will now reformulate 
the problem employing the language of Berry phases. Namely we will adopt a geometric ap- 
proach to evaluating the crystal momenta of magnetic systems ll 1811 . and show that the momentum 
of eq.©, which was central to the LSM argument re-emerges as an additional Berry phase which 
is generated upon application of the twist. An incarnation of this geometric effect will arise in the 
form of vortex Berry phases in the low energy field theory described in section [TV] which can be 
constructed for arbitrary d. 



B. Berry phase argument 

We begin by quickly recalling how a spin Berry phase typically arises. Consider a spin-S* object 



represented by a spin coherent state[28] |n(t)), where the unit vector n(t) specifies (in a semiclas- 
sical sense) the spin's orientation. When the dynamics of the system is such that n(t) undergoes 
an adiabatic rotation, returning to its initial value at the end of the excursion, the wavefunction 
accumulates a net phase of S § dn(t) ■ a(n(t)) = Su>[n(t)]. Here a(n(£)) is a Berry connection 
defined by (n(i)|n(£) + 5n) = exp [iS&(n(t)) ■ Sn), C is the loop on the unit sphere mapped out 
by the trajectory {n(i)}, and cj[n(i)] the solid angle enclosed by C. 

One can conceive of situations where a similar phase accumulation is induced along a spatial 
(as opposed to temporal) extent of a many-spin system by a gradual spatial change of n. Such an 
example was elaborated by HaldaneL18] in his study of the crystal momenta of a ferromagnetic 
spin chain. Here one is concerned with an instantaneous spin configuration written as a direct 
product of spin coherent states, 

N 

|{n,-}) = (g) In,} . (4) 
i=i 

As in the previous subsection, information on the crystal momentum P is gained by inspecting 
how the generator of translation T = e lPa (a: lattice constant) affects this configuration. However, 
we should keep in mind that the state © is in general not an eigenstate of T, and we will in this 
semiclassical approach be evaluating the expectation value ({n.,}|T|{n.,}) instead of the eigen- 



value itself fl 1 8h . Since by definition ({njHTKn,}) = [7 = 1 (n.j|n :? -_ 1 ), (we assume a periodic 
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boundary condition (PBC), where the (iV+l)th site is identified with site 1) we find, in analogy to 
the familiar temporal Berry phase, that in the continuum limit 

({n(x)}|T|{n(x)}> = e iS ^ n ^. (5) 

The solid angle cj[n(x)] is now associated with the loop traced out by the snapshot configuration 
{n(x)}. Eq.© suggests that the crystal momenta of a d=\ ferromagnet is a topological quantity 
which is best described as a quantal phase defined modulo 2n. 

As illustrated in the previous subsection, evaluating the crystal momenta carried by the low- 



energy states is the central step in a LSM-type scheme[15]. It is therefore tempting to derive the 
counterpart to eq.© for our present problem, since by introducing such a relation into the LSM 
program, we can expect to arrive at a geometrical picture underlying the magnetization properties 
of antiferromagnets. We show below that this is indeed possible. In doing so however, the follow 
points require clarification. (1) The ground state of a quantum antiferromagnet is generally com- 
plicated, and cannot be expressed in the form of a spin coherent state as in eq.@. However, one 



can show Il29h that if we take the large-S 1 limit, the quantum ground state will generally approach 
a conventional collinear Neel state. This is the basis of our use of a coherent state ansatz in the 
following semiclassical treatment. (We will explore in appendix [A] what will happen when we 
relax this ansatz; see the remark of the end of this section.) (2) It is also pertinent to observe that 
eq.© is sensible only when the spin orientation n(x) is of a smoothly varying nature. Meanwhile 
antiferromagnets in an external field would generally have components which vary on the lattice 
scale. For the Berry phase framework to work, therefore, we must devise a way to represent our 
system in terms of a coherent state labeled by some slowly-varying field. 
Let us start then from a canted configuration 

iij = ((— iy cos <j)j sin 9j, (—iy sin <pj sin 9j, cos 9j), (6) 

where the unit vectors {n,} represent the orientation of spins partially polarized by a magnetic 
field applied along the z-axis(see FIG. [I])- We can identify the unit vector 

Nj = (cos (f>j sin 9j , sin <fij sin 9j , cos 9j ) (7) 

as a slowly-varying (unstaggered) vector field for which the corresponding Berry connection 

Sa(N(x)) = ({N(z)}|HV N(s) )|{N(x)}) (8) 
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FIG. 1: Schematic view of the slowly varying angular variables {4>j}, and the magnetization per site m. 



is a valid construction. This motivates us to introduce the unitary operator U = exp 



which transforms ({nj) into \{Nj}), i.e. U\{nj}) = \{Nj}). Noting that fuf- 1 



exp 



U (with PBC and iV=even assumed), we find that 

(K-}|T |{ % }> = ({n^le-^-^-^^^fC/lin,-}) 



~ e i7r Sf=i( s - 5cos ^)({N i }|f|{N,-}}, (9) 



where the last line is a large-S* result and is a consequence of the minimal uncertainty property of 



the spin coherent state yOD. We may now incorporate Haldane's formula, eq.© to obtain 

({nj}\f IK-}) = e ™?$=i(S-S<x*O s ) e iSa,M*)]. (10) 

Thus we have separated out the intrinsic part, where the Berry phase appears, from a momentum 
offset arising from the staggered nature of the ground state. (We pause to observe that in the 
absence of the magnetic field, a similar attempt to extract a Berry phase associated with smoothly 
varying components would run into difficulties, since that would involve, instead of a simple 
unitary transformation the simultaneous flipping of all three components of spins residing on every 
other site, i.e. time reversal operations.) The existence of this magnetization-dependent offset is 
also easily seen from the exact Bethe-Ansatz solution! 3JJ] of the S = 1/2 XXZ chain. 

The next task is to compare eq.dTOb with the corresponding expectation value for a LSM-twisted 
state. The latter state is defined by |{n^ SM }) = Uisu\{^-j}), where U hSM = e ^=^^ . We 
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note that [U, ?7lsm] = 0. (Following OYA we are assuming that the interactions are sufficiently 
local and the Hamiltonian possesses a rotational symmetry around the z-axis.) The LSM-twisted 
counterpart of eq.dlOl) is 

({n^ SM }|f |{n^ SM }) = e ^T,f=i(s-Scose ] ) e iSuj[N^ M (x)]^ ^ 

where N LSM (a;) differs from N(x) by a shift of the azimuthal angle, i.e. <Plsm(x) = <f>(x) — ^-x, 
L = Na. Notice that the twist leaves the offset portion of the momentum unaffected. The right 
hand side expressions of eqs.dTOb and (fTTI) will coincide if 5'a;[N LSM (x)] = 5a;[N(x)] + 2im, 
where n G Z. Following the usual logic of LSM-type arguments we may interpret this as the 
semiclassical expression of the necessary condition for the occurrence of a spin gap, i.e. a mag- 
netization plateau. We can make contact with the OYA theory by using the spherical coordinate 
representation for the solid angle, 

Scj[N lsm (x)] = S I dx{l - cos9{x))d x (j) LSM {x) 
Jo 

27T S 

= Slu[N(x)} — / dx(l-cos9(x)) 



= Su[N{x)]-2ir{S -m), (12) 

where m = | f» dx cos 9 is the magnetization density. This is the Berry-phase derivation of the 
LSM momentum shift ©. The aforementioned condition thus translates into the quantization 
rule S — m 6 Z. This argument is readily extended to the case where the unit cell consists of 
r > 1 sites; the relevant quantity then will be the expectation value of T r , leading to the spin-gap 
condition rS , co'[N(x)] = r5'u;[N LSM (a;)] + 27m, n e Z, which in turn yields the OYA quantization 



ruleBalrJ r(S-m) e Z. 

Finally a word on gauge independence is in order, as eq. (IT2|) involves a particular gauge choice 
for the monopole vector potential a(N(x)) (the Dirac string goes through the south pole). One 
finds that relocating the string to the north pole merely shifts the crystal momentum by AnS, 
which is immaterial. Likewise, other gauge choices consistent with the spherical geometry of the 
target manifold will leave the results unaltered. 

As noted above, we have assumed on semiclassical grounds that it suffices to deal with ground 
states (and their twisted counterparts) which can be expressed as a spin coherent state. One may 
wonder though whether we can extend the present argument to more generic ground states which 
are superpositions of coherent states. We show that this is indeed possible in appendix lAl 



IV. EFFECTIVE ACTION AND DUAL VORTEX THEORY 



The findings of the previous section tell us that for d = 1, a simple semiclassical picture 
involving partially polarized spin moments is capable of accounting for the results of OYA, once 
we realize that certain quantum interference effects are at work. The latter are naturally described 
in terms of Berry phase concepts. 

This prompts us to seek a generalization of such geometric interpretations to arbitrary d. Since 
the LSM approach is now unavailable, we will shift gears and carefully work out an effective low 
energy theory valid in arbitrary d, invoking once more the semiclassical picture of canted spins. 
We will find that Berry phase effects are again present, now manifesting themselves in the form 
of quantum interference among vortex configurations. These interferences can drastically alter the 
low energy properties of the system. (We will come full circle by later specializing to d = 1, and 
identifying a common root that this effect shares with the LSM-type argument presented in section 
Hill ) We also discuss possible ways to perturb the system so that the Berry phases (and hence the 
magnetization properties) are modified. 

A. Effective Action 

We now substantiate the foresaid by deriving a low energy effective action for an antiferromag- 
net coupled to an external magnetic field. The effective theory, summarized in eq.(l23l) below, is 
valid irrespective of the dimensionality. 

For the sake of concreteness, we will consider the following Hamiltonian: 



on a d-dimensional hypercubic lattice. The first term corresponds to the usual Heisenberg Hamilto- 
nian and the second to the single-ion anisotropy. The external field H is applied in the ^-direction. 
If we regard the phase of S + and the local magnetization S z respectively as the Josephson angle <p 
and the boson density n, it is easy to see the similarity to the Josephson junction array: 



or, in more general terms, the boson Hubbard (BH) model. In the large-S 1 limit of the model 
(TT3T) . spins assume a canted (or,conical) configuration with the XY components aligned in an anti- 




(13) 




(14) 
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parallel (antiferromagnetic) manner: 



( Js 2 -m 2 j cos(Qtj) \ 



S 2 — m 2 sin(Q-rj 



V 



3 

m 



(15) 



/ 



where Q = (n/a, . . . ,it /a). When either the field H or the single-ion anisotropy D is finite, a gap 
opens at k = Q and the only gapless mode (a transverse spin wave) appears at k = 0. 

Below we extract a low-en ergy ef fective theory from the model of eq.(fT3l). incorporating stan- 



dard path-integral procedures ||32l. 



33|]. We illustrate this for the one-dimensional case. (Generaliz- 



ing to higher dimensions is straightforward.) Guided by the classical solution, we parametrize the 
fluctuation around the above canted state as: 

n ±( x ) = (_i)r e ±W*) ( S in0(°) + 59(x) cos# (0) ) 
n z {x) = cos(# (0) + S9(x)) « cos# (0) - 89{x) sin0(°> , 



(16) 



where 9^ = cos^ 1 ( js^j+d) ) " ^ ne two ec l uat i ons above each correspond to modifying e* Q 
in ((131) as ^Q-rji^M^ an( j expanding rnj around the average m = S cos9 ( -°\ An inspection of the 
classical Poisson bracket relation for the spin variables suggests that we identify the following as 
a pair of canonical variables (a denotes the lattice constant): 



q(x) = (p(x) , p(x) = —S sin9^°'59(x) = all( 



x) 



which satisfy the equal-time commutator[34] 



[(p{x), Tl(x')] =i5(x-x'). 



From the expression 



S„ ~ S cos 9 (0) + a U(x) + ■ ■ ■ = m + aU(x) + ■■■ 



(17) 



(18) 



(19) 



it is obvious that II describes the longitudinal fluctuations around the average magnetization m. 

Casting these into path-integral form and retaining terms up to second order in cp and II we 
obtain the action S c i + Sbp where: 

H 2 



S d = jdrjdxa(2J + D)^(r)- 1 -^(l ^ {D + 2J) 
Sbp = iS (1 — cos 

dxdr d T tp — i dr dx d T (pH 



a j dr J dx (d x ip)" 



(20) 
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The two contributations <S c i and Sbp each come from the classical Hamiltonian and the sum over 
the Berry phases of each spin. Although the first term of Sbp is a total derivative, it cannot be 
dropped from the effective action for a reason which will become clear below. Following this 
intermediate step, we integrate out the Il-field to obtain the effective action for the angular field cp: 

S eS = Jdxdr^(d T ^ 2 + ^ad x ^ 2 ]+t^^ Jdxdrd T ip, (21) 

where the spin stiffness ( and the spinwave velocity v are given by 



S V A(D + 2J) 2 S 2 )' Y 2J(D + 2J) K ' 

Additional interactions will merely renormalize the values of £ and v. Note that the 92-field appear- 
ing in (l2"TT) is defined on a universal covering space of a circle and the last term counts the winding 
number of the space-time history. In section |V] we will show that it is also possible to arrive at 
this term by carefully summing over the spin Berry phase terms iSuj[n(ij)] associated with each 
lattice site. 

In the above, we have assumed that in the transverse (XY, here) direction, the slowly-varying 
degree of freedom is a staggered component with wavevector Q, i.e. possesses at least a short- 
range antiferromagnetic order. However, the argument goes exactly in the same manner for the 
spiral (helical) magnets as well; the external field kills two of the three Goldstone modes and the 
remaining one is described again by </?. Therefore, our effective action (|2T|) is applicable equally 
well to non-frustrated- and frustrated cases. 

Our derivation of <S e g ((2TI) for d = 1 readily generalizes to any spatial dimension d, and leads 
to the following effective action: 

S eS [(f(T, x)] = «S t op[<p(r, r)] + Sxx[<p(r, r)], 



Stop = i / drd rpd T (p 



Sxy = / drd d v 



^(^) 2 + ^(Vy,) 2 



(23) 



where p = ^Jr, K T = (/v 2 and K± = (. Eq.(|23l bears the form of an XY model in (d + 1) 
dimensions supplemented with a topological term. An identical action was previously employed to 



describe the hydrodynamical properties of superfluids [|22L 1231. 124|] . In this analogy, p plays the role 
of a uniform offset value of the superfluid density, as one can read off from eq. (IB3l) of appendix[B] 
There is also an apparent similarity to the low energy theory for the BH model at incommensurate 
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filling factors H2p|. 2111 . which is natural in view of the correspondence between eqs.CP ) and (14); 
we will come back to this connection later. As already mentioned, the topological term S top is a 
total derivative and as such will not affect the classical equation of motion. It does however have 
profound influences on low energy physics when the quantum effects of space-time vortices, to 
which we now turn, are properly accounted for. 



B. Dual Vortex Theory 

The sensitivity of the topological term to vortices is already apparent from the simple observa- 
tion that a phase winding (in imaginary time) of 2% at spatial position rj will yield a nonvanishing 
contribution 5S t0 p = i2np. To understand the consequence of such effects, it proves convenient to 



apply to eq.(l23l a standard boson-vortex duality transformation ll35ll and recast the action in terms 
of vortex events: point-like space-time vortices (phase-slips) in (l+l)d, vortex loops in (2+l)d 
(world-lines of point-like vortex- anti vortex pairs), and closed vortex surfaces in (3+l)d (world- 

ceeping 



sheets of vortex rings). For brevity we will generally resort to continuum notations while 
track of the lattice origin of our model. Following steps well- accounted for in the literature A3 5L 13611. 
we extract the following action (derivations are supplied in appendix IB]), 

^vortex = ^Coulomb 4~ <5gp , (24) 

where iS>c ou iomb re P resents the (d+l) dimensional Coulombic interaction among space-time vortex 
objects mediated by the kernel 

;A) 2 + 77^(V) 2 1 , (25) 



-d 2 J \2K ± y ' 2K T 

and iS>Bp tcx is the Berry phase associated with vortex events. The latter inherits the information 
on the "superfluid density" p oc S — m contained in the original topological term and takes the 
following forms: 



cvortex 
°BP 



2it(s— m) \—\ (0) /-i , in j 

-V^EiW ( 1+1 ) d 



"-^E^Al ( 2+1 ) d (26a) 

(o) 
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where is a solution to 



'_ e da i0) = 6 


(l + l)d 




< e d a (0) = 5 


(2+l)d 


(26b) 


-e 8 a {0) -8 


(3+l)d . 





By using the explicit solution to eq. (|26bl) . S^ tex can be written compactly as: 

;^^£,^ (i+i)d 



ovortex 
°BP 



^■T.M? (2+i)d 



(26c) 



The summations on the right hand side are to be taken over all vortex events, and the q/s are their 
vorticities. Xj denotes the spatial coordinate of the j-th (l+l)d vortex, A xy the area bounded by 
the projection onto the xy-plane of the j-th vortex loop Cj in (2+l)d. V* yz is the volume of the 
j-th vortex surface in (3+l)d projected onto the xyz space (i.e. the net real-space volume occupied 
by a vortex-ring through its lifetime). 

The implications of the vortex Berry phase factors are clear. In (l+l)d, the profile of vor- 
tex events, when projected onto the spatial coordinate axis consists of points residing on dual 
sites. The spatial coordinate difference of any pair of vortex events (not necessarily occurring 
at equal times), Xj — Xj, are therefore integer multiples of a. Likewise, in (2+l)d and (3+l)d, 
A xy = a 2 x integer and V xyz = a 3 x integer, respectively. Hence when S — m Z, vortices will 
create contributions to the partition function that are weighted by oscillatory phase factors, gen- 
erally leading to a destructive interference (unless events with vorticities of equal magnitudes and 
opposite signs are confined in pairs). (For the moment we are leaving aside the commensurabil- 
ity effects which set in when S — m is a rational non-integer number.) In other words vortices 
are, under this condition unable to contribute to the partition function. In (l+l)d, for instance, a 
Kosterlitz-Thouless transition into the plasma phase is prohibited quantum mechanically regard- 
less of the value of the coupling constants. 1 32. 3711 When S-m G Z, in contrast, the phase factors 
trivialize and the partition function reduces to that of the (d + 1) dimensional XY-model. Vortices 
are then able to proliferate and condense, if energetically favorable, and drive the system into a 
disordered state. 

With the hindsight of the results just mentioned, we briefly reflect upon the LSM approach of 
the previous section and attach to it a simple physical picture. Given that the system in question can 
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FIG. 2: A vortex loop event Cj in (2+l)d. Projected onto the xy plane, this corresponds to a vortex- 
antivortex pair created at point 1, which pair-annihilates at point 2. 

basically be regarded as a (l+l)d XY model subject to a periodic spatial boundary condition, we 
shall consider a one dimensional superfluid (with superfluid density p s ) in a ring geometry, a setup 
often employed to discuss the stability of persistent supercurrents. The LSM twist in this language 
is nothing but a phase slip process which causes the total phase difference along the ring to change 
by ±27r. It is well known[32, 38] that such an event will generate a Galilei boost, i.e. an increment 
of the center-of-mass momentum of the superfluid by the amount AP = J dxpd x (p = ±27ip. 
The latter leads to a full agreement with the findings of section |Tll] provided that we are allowed 
to equate p s with the coefficient of the topological term, p = . As mentioned earlier, a short 
inspection of eq. (|B3l) in appendix|B]confirms that the topological term indeed functions as a source 
term for the superfluid density, imposing precisely this value. In this way we find that the vortex 
Berry phases displayed in eq. (!26cl) and the momentum increment due to the LSM twist share a 
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common origin [|39h : the topological information encoded in the term «S top of eq.(l23l). Within the 
present physical picture, the novelty of the condition S — m E 7h may be understood by observing 
that it allows for the extra phase winding due to the phase slip to occur without the need for the 
superfluid to transfer excess momenta to phonons. (If not for this condition, such transfer will be 
generically forbidden at zero temperature.) In other words, for this special case it is possible (if 
energetically favored) for the ground state to be a superposition of different circulation numbers 
t^t § dxd x (p G Z, i.e. a condensate of phase slips. 

While we have so far treated the vortices at a "first-quantized" level, the Berry phase factors 
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4JJ], which 



of eq. (|26cl) are readily incorporated into second-quantized vortex-field theories 112 ll. 
provide a framework more amenable to detailed analysis. Referring to the later half of appendix 
|B]for details, we display here their main structures: 

$^[01 = Skin[$ +V J d 2 r cos[27r(S - m)x + 2vr$ (27a) 
SlTK X] = S^nK) + E ^ / * r C ° S IM\X -\ ~ MS " m)x)} (27b) 
5 d 3 u + ai )d [^' c m] = <WM + E J d A r cos[2vr(A At c 1/ - A„c M ) - 2^(6^ + S^S^S - m)x)\ 



(27c) 



The kinetic terms denoted by 5 kin are quadratic; S^ 1 ^ is a gaussian theory, while the higher- 
dimensional counterparts are Maxwellian terms associated with the dual gauge 1- and 2-forms, 
and b^ u . In principle, the action can be supplemented with additional terms consistent with the 
symmetry requirements of the problem. As explained in appendix El d x {p (d=l) and the gauge 
curvatures e T ^ u d^b u (d=2), e Tflu xd^b u x (d=3) are directly related to the local magnetization (or in 
the superfluid analogy, to the superfluid density). The fields x an d c M are introduced into the 



theory to implement the continuity of the vortex current (2+l)d/vortex-loop current (3+l)dll21 
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4lh . For each case the Berry phase effect enters as a spatial modulation 2n(S — m)x within the 



cosine term, again causing oscillatory behavior unless the forementioned condition is met. Similar 
expressions (including higher harmonics terms) have been discussed in (l+l)d using bosonization 



methods. P, 



17] 



It is worth noting that for (2+l)d, the above expression is identical in form to the free en- 
ergy of a (classical) 3d lattice superconductor, subjected to an external magnetic field of strength 
2tt(S — m) (expressed in the Landau gauge). The problem at hand therefore reduces in this case 
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to the energetics of a system of Abrikosov vortices immersed in a dual lattice superconductor; if 
these vortices manage to destruct superfluidity in this dual theory, it is tantamount to a massless 
superfiuid state of the original XY model. If on the other hand superfluidity remains intact in 
the dual theory, the original theory is disordered, i.e., is massive. This observation enables us to 
make contact with a related argument due to Lee and Shankar[42] (LS), who also employ a duality 
technique (which is somewhat different from ours, and in particular do not explicitly emphasize 
the relevance of Berry phases) to map quantum 2d lattice systems into an ensemble of Abrikosov 
vortices. Following LS, it is natural to expect that when S — m is an integer (or more generally, 
is a rational number-a case we had not explicitly addressed up to now) the ground state is a peri- 
odic vortex lattice commensurate with and pinned by the underlying lattice structure of our spin 
system. For irrational S — m, they will form a floating lattice and destroy the superfluidity of the 
dual superconductor. 

The case where S — m is a non-integer rational number, should however be treated with special 
caution, since subtle quantum effects may lead to novel physics which the naive superconductor 
energetics can miss|l9|]. To illustrate what can happen, we go back to the vortex Berry phases 
of eq. (!26cl) . and put S — m = p/q, q > 2. Singly quantized vortices/antivortices (qj = ±1) are 
clearly frustrated; they will suffer the accumulation of an AB-like phase of ^ upon encircling a 
dual plaquette. This leads to the destructive interference that we had been describing earlier in this 
section. Vortices whose winding numbers are integer multiples of q, on the other hand are free of 
such effects. This leads us to expect that a highly unusual phase can appear, in which only vor- 
tices with vorticity ±q are able to condense; such a tendency can be stabilized by the q-th higher 
harmonics of the cosine term of eq. (|27l) . which is the leading non-oscillating term in the harmonic 
expansion. (Needless to say this commensurability effect becomes weak when q is too large.) 



Since a condensate of higher- winding number vortices can sustain fractional excitations [43 



this simple argument points to the possibility of yielding a fractionalized phase in antiferromag- 
nets in an applied magnetic field. Actually an even richer variety of phases are possible in this 
commensurate case, and can be investigated in its full generalit y by adopting methods initiated by 



26, 



46, 



470 . It is readily seen that 



Lannert et al |45|], and further developed by other authors [|25L 
there are actually q degenerate low energy modes within the magnetic Brillouin zone for frustrated 
vortices hopping on a square (dual) lattice under the influence of a —-flux piercing each (dual) 
plaquette. The natural procedure would thus be to construct a Ginzburg-Landau-like action involv- 
ing q species of vortex fields, and to study the condensation of various composites of such objects. 
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We will not work out the details of this approach here as they can be found in the literature[25] 
(albeit in different contexts). In section[V]we will instead focus on the important case where S — m 
is half-integer-valued, and illustrate the nature of possible fractionalized phases which can emerge. 



C. Correspondence with the Boson Hubbard Model 

It is instructive to understand how the above conforms with properties of the well-known chem- 
ical potential vs. tunneling parameter phase diagram of the BH model[4] . The most generic feature 
of the BH phase diagram are lobes of Mott insulator phases lined up along the chemical poten- 
tial (/i) axis in the weak tunneling region. The state within each lobe, being incompressible, is 
characterized by a constant integer- valued boson density N. As evident from the foregoing (see 
also the discussion of SachdevyJ]), there is a set of correspondences here with the present prob- 
lem which may be summarized as: fi <->• H, N <-> S — m, Mott insulator^ spin gapped phase, 
superfiuid phase^gapless phase. The BH counterpart to our problem of determining the gapful- 
ness/gaplessness of the antiferromagnet at fixed S — m thus amounts to examining, at a fixed value 
of N, whether it is possibility for a superfiuid to directly enter a Mott insulator phase upon the ap- 
proach to the weak tunneling regime. Here we would like to place this analogy on firmer grounds 
by recalling how Berry phases affect the critical properties of the BH model. 

It is known[4] that when N £ Z the constant- N contour in the superfiuid phase merges with the 
Mott insulator phase (with the same N- value) at the tip of the lobe. Meanwhile the contours for 
noninteger-iV will disappear into the region in between two adjacent Mott insulator lobes; the fate 
of these contours (especially those representing commensurate fillings) is a delicate matter[48] and 
will ultimately depend on the short-range physics. The former situation, in which a superfiuid-to- 
insulator transition is generically possible, is therefore reminiscent of the S-m£Z case of the 
antiferromagnets. The reason behind this similarity can be traced to the behavior of the topological 
term present in the low energy action; namely, both systems correspond to cases where this term 
ceases to be effective. We have already discussed in length the irrelevance of Berry phase effects 
for the S — m E Z antiferromagnets. Meanwhile the topological term for the effective action of 
the BH model is strongly constrained by a U(l) gauge symmetry which the original Hamiltonian 



possessedlll, |21fl. This constraint can be used to show that the coefficient of this term vanishes 
identically right at the lobe tip. Hence the quantum phase transition which takes place along the 
integer- N contour at this point is in the universality of the (d + 1) dimensional XY model, which 
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is the same universality class relevant to the antiferromagnets with S — m G Z. 
D. Effects of antisymmetric interactions 

In closing this section we briefly discuss a possible way of modifying the vortex Berry phases 
we have discussed above. Consider, in the superfluid analogy, allowing a bulk supercurrent J 
to flow through a rf-dimensional sample; this corresponds to perturbing the spin system with a 
Dzyaloshinskii-Moriya (DM) type interaction 

d 

n DM = 5^ A^(S r xS r+ eJ . (28) 

r a=l 

It is easy to see that this perturbation induces a non-zero spin current (which translates into a 
supercurrent), and can be accounted for in the effective theory of eq. ((2~3~l) by simply adding the 
term iJ a V a y2, where J a ~ D a is the induced supercurrent. This addition will obviously effect 
the vortex Berry phases. For instance, while the (2+l)d expression in eq. (|26c| ) basically counts 
the number of bosons residing within the shaded region of Fig. 2, we must now correct for the 
migration due to the flow. (Physically, this will cause a change in the quantum dynamics of each 
vortex. Besides the Magnus force (perpendicular to the velocity of the vortex) coming from the 
original topological term, vortices will suffer an additional "Lorentz force" which is perpendicular 
to the supercurrent.) This shifts the Berry phase term by 

A^ P +1)d = *(2vr) 9j(J,Af + J y A] x ) , (29) 

j £ loops 

where A^ T and A™ are each the area of the projected image of Cj onto the yr- and rx planes. A 
similar consideration leads to the following Berry -phase shift in (l+l)-dimensions: 

AS$£ 1 * = i(2ir)J x '%2q j T j . (30) 

j 

Now let us consider how the non-zero supercurrent J a (or, non-zero DM interactions) affects 
the magnetization process for the simplest (l+l)-dimensional case. In accordance with eq.(l27l). the 
low-energy physics in the vicinity of S — m G Z may be described by the following sine-Gordon 
model: 

< 5 duai )d [$ = <W^] + y J d 2 r cos[2vr(5mx + J x r) + 2tt^] (31) 
where 5m measures a small deviation from S — m G Z. 
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By a straightforward extension of the methods used in Refs.37|, |49L we can carry out a 
renormalization-group analysis. According to the 1-loop renormalization-group equations, there 
are two different length scales: (i) a length scale £(<5m, J x ) = 1/yJ (5m) 2 + (J x ) 2 set both by 5m 
and by J x (or, D) and (ii) another set by the lowest particle (i.e. 5S Z ^ 0) excitation gap A so i 
in the absence of 5m and J x (if it is finite). When A so i = 0, nothing competes with 5m or J x 
and the system is in the gapless (Tomonaga-Luttinger liquid) phase. When A sol ^ 0, on the other 
hand, we have two different possibilities depending on the ratio between £,(5m, J x ) and A~J; if 
£,(5m, J x ) < A^p renormalization stops well before the system starts feeling the effect of the 
pinning potential and the system flows into the same gapless phase as above. If £,(5m, J x ) > A^J, 
the system is renormalized into another phase. The nature of this phase may be analyzed by map- 
ping the system to a system of (nearly) free fermions[50]; we find that the gap responsible for the 
plateau is robust against J x , i.e. the DM interaction D. 

V. NATURE OF PLATEAU PHASES - Z 2 GAUGE THEORY 

In the previous section, we had seen that for S—m e Z vortices do not suffer from interfer- 
ence effects due to the Berry phase (I26cl) . As a consequence, the usual (d+l) -dimensional XY 



transitionll4i 15 It! separates the gapless XY-ordered (superfluid) phase at \arge-K XT from the small- 
er plateau (i.e. insulating) phase where vortices proliferate. For generic incommensurate values of 
S — m, on the other hand, destructive interference among different configurations of topological 
(vortex) excitations should lead to quite different small-i^ behaviors. However, it is not easy to 
capture the fine structures which may emerge as we approach the limit K — > 0. For this reason, 
we utilize a slightly different formulation of the problem[26] more suited to study the strong cou- 
pling phases, and attempt to see more closely how the topological excitations destroy the gapless 
superfluid phases. 

The derivation presented below, leading to our effective theories rely heavily on methods de- 
scribed in ref 26. We are lead however to interesting differences which we will highlight as they 
appear. Let us discretize the (d+l) -dimensional space time and put the quantum XY model (|23|) . 
using the usual Villain form, on the (d+l)-dimensional hypercubic lattice: 

d+l 

Sxy = -K XY - 2nm ^) 2 ■ ( 32 ) 

The XY coupling K X y is proportional to K Tj ± in eq.(l23l). The above XY action should be supple- 
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mented by the Berry phase term: 



S S p = i(2S)^2A jtT , (33) 

3 



where Aj !fl (fx = 1, ... , d+1; x^+i = r), the lattice version of the spin gauge field (or equivalently 
the CP 1 gauge field), is half the area of a spherical triangle enclosed by a triad of unit vectors 
n (a fixed reference vector which can be chosen arbitrarily), n(j) and n(j + ft). Unlike the 



5, 



BP 



corresponding expression for 5 = 1/2 antiferromagnets in the absence of magnetic field 
here does not contain an sign- alternating factor, reflecting the canted nature of the classical ground 
state. In addition to these, we add (by hand) the simplest 'Maxwellian term' allowed both by the 
arbitrariness of n and by the 2tt -redundancy (compactness) of Aj >fl . The (l+l)d expression (we 
will be dealing with this case in most of the equations to follow), for example, is 

^Maxwell = i E (V A /A> ~ 27 %* f ■ (34) 

In the above equation, the integer field qj* resides on the dual lattice site. Carrying out the Poisson 
resummation and Gaussian integration, we obtain: 

3* 3 

It is convenient to rewrite the Berry phase term, eq.(T33l, as 

S B p = i (2S) V E A> 



3 A* 



\ (0) 



(36) 



= i (25) E( 6 ^ A A"4?>; u '7* 

where we have introduced an off-set gauge field which satisfies e Ml/ A z/ a^*' ) = 5^ T . This can be 
solved explicitly as, e.g., 

af) = -j* x . (37) 



As an important physical aside, we note that for the canted spin configuration of eq.(fT5l). the 
lattice curl e^A^Aj^ (which usually represents the antiferromagnetic spin chirality fluctuations) 
is proportional to the vorticity e^A^m^ of the XY spins on a plaquette surrounding j*. (This is 
similar to the situation of the easy plane antiferromagnet in the absence of a magnetic field treated 



in ref. 



261) In terms of the latter, the Berry phase of eq.(l36l) reads: 

5 BP = i 2ir(S - m) ^(e^ u A^m jjV )afj . (38) 
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Notice that if we identify e^A^rrij^ with qj of the previous section, the above Sep exactly coin- 
cides with the first line of eq. (|26cl) . This is readily generalized to higher-dimensional cases. For 
example, in (2+l)d we may write 



W 0) 



(39) 



.14' 



with (afj x ,afj y ,afj r ) 
result (I26cl) . 



-j*/2,j*/2, 0). It is easy to verify that this reproduces the previous 



Rescaling a,-* i— > (2S)aj* in ^M^xweii an ^ collecting the terms S X y, ^M^weii an ^ <^bp, we arrive 
at the following action: 

p2n r i 

Z (i + i )d= Il^ ex P -2^Y^(A^,-2vrm,, M ) 2 



<7n 



(40) 



3* j 

To make further progress we now need to fix the coefficient of the Berry phase term. Here a 
crucial difference arises with eq.(53) of refj^, which corresponds to setting m = (and S = 1/2) 
in our eq.(l40l). By varying m, we are able to probe through a variety of different effective theories, 
each characterized by a different set of vortex Berry phase factors. Below we will inspect a few 
representative cases. 

When S — m £ Z, the U(l) gauge field aj* and the vortices decouple from each other. The 
a-summation in eq. (|40l) can then be carried out trivially, to yield the classical XY model in (d + 1) 
(d-l here) dimensions: 

K XY 



e j n^ exp E( a ^ " 2j[m ^ 

m,,u} ° j L 3,1* 



{ m j>} ' 



(41) 



Thus in agreement with our arguments of section HVB I we expect to encounter the usual classical 
XY transition triggered by vortex proliferation as the XY-coupling K X y is varied. 

A more interesting situation arises for S" — m G Z + 1/2. Here the partition function (1401) 
depends on the parity of the link variable m JiA1 , which leads us to introduce the Ising variable 



through||26ll 



m 



1 



(/?> e z, Sjj+p, — ±i) . 



(42) 
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Plugging this into eq.(l40l) and carrying out the summation over the a 3 * explicitly, we obtain[52] 
the following partition function: 

/2tt 

j J dtp j q- s ^- s ™«™ , (43) 

where the action of the quantum Z 2 gauge theory 5i s i ng and its coupling to the matter (XY) field 
each take the form: 

5 Ising = -K 1Gr (g m ) ( II s i-i+A ) + ™ E ' ? ? ( 44a ) 



□ \ n 



^matter = -4-^XY E S W+A COS Q^M^iJ ' 



(44b) 



J4> 

The effect of the Berry phase term of (l40l) is now encoded in the second term of iS>i s i ng . The gauge 
coupling Ki G j, a monotonically decreasing function of g, is defined through the relation: 



e 2K 1GT (g) 



E e ^ n VE( _1 ) v ^ 2 ' (44c) 



The appearance of the factor 1/2 in the cosine of (I44bl) is worth noting; it indicates that the 
Ising gauge field couples to fractionalized bosons, whose creation operators can be identified as 
fet ~ . A theory with the same basic structure, i.e. fractionalized bosons held together by 
Ising gauge fields, is obtained in (2+1) dimensions as well. Taking the limit .Kigt=0, an explicit 
summation over the gauge variables Sj t j + p, leaves us with a theory in which coefficients are dou- 
bled wherever (A^ipj) arises; hence we can expect the fractional bosons to be confined in this 
limit. A more nontrivial phase may emerge by looking into regions with larger Ki G T-values. It 
turns out in fact, that deconfinement can occur in (2+1) dimensions for large -Kigt* realizing a 
fractionalized plateau phase with topological order. Such details are most conveniently analyzed 
by going back to the form of eq.(l40l) and performing, in similarity to the previous section, a duality 
transformation. 

Apart from the difference in the expression for af), the derivation of the dual theory proceeds 
according to the prescription of ref]26l The final form for (l+l)d reads: 



{a*}{n*} L j* n,v 

(45) 

where {nj*}, whose lattice curl is the boson density, is an integer- valued field defined on the dual 
sites. The kernel M?'?, governing the density-density interaction reduces to M^' v , = SjfS^ for 
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the simplest case of spin systems with short-ranged interactions; here it is written in a form which 
accounts for a more generic Hamiltonian[53]. A similar dual action is also obtained for the (2+l)d 
case. While we have displayed the form of the dual theory for general S — m, we continue to 
concentrate here on the case 5-m6Z + l/2. 

(a) & — 4 J- — 4 4 — 4 4= 

(b) ^ J^_^ 

(c) a a 4 ^ ^ 4* 

kink 



FIG. 3: A typical spin state found at an m = 1/2 plateau in spin-1 chains (S — m = 1/2). Lattice translation 
symmetry is broken spontaneously and as a consequence we have two degenerate ground states (a) and (b). 
Half (arrows) of each spin-1 degree of freedom is frozen by a strong magnetic field and the remaining 



fluctuating half forms valence-bond solid states[55]. The lowest excitation is a kink which carries an S z 
quantum number of 1/2. 




FIG. 4: Two-dimensional analogue of the spin-1 partially-polarized valence-bond-solid state shown in 
FigfJ] Again, the unpolarized 'fractions' of the spin-l's form a valence-bond pattern. Dashed lines denote 
domain walls separating four different valence-bond patterns. At the core of the Z^-vortex is a fraction- 
alized excitation carrying S z = 1/2 (large arrow). The resemblance of the spatial patterns as well as the 



fractionalized nature of the core excitation with those of Refs 



56. 



57 is evident. 



Insight into the phase structure of Z ld (and Z 2d ) comes from following the analysis carried out 
in refko. We begin by considering the limit g m — » oo. Under this condition, the fluctuations in 
the a,j* -field are quenched aj* > a^J and the model (|45T) reduces (for a diagonal M^) to the 
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dual vortex model of appendix |B] (compare with eqs. (IB3l) and (1B4I) ). (Conversely, eq.(l45l) extends 
the standard dual- vortex model treated in appendix |B] to allow for fluctuations of the background 
boson density.) The latter is expected to be in a superfluid (or, XY) phase I148n . which persists into 
the region K X y 3> 1 • 

Turning next to generic values of g m , the aj* -field starts to retain its dynamics and induces fluc- 
tuations in the background boson density (spin density in the present physical context) e T \A\dj* . 
These fluctuations may stabilize a new phase by forming bond-centered orders. (By retracing our 
steps backwards through the sequences described above, one sees that the g m -terms of eq.(l45l) and 
its (2+l)d counterpart are associated with the energies of the spatial links of the direct lattice). In 
fact, by recasting our dual theory into a 2d height model[26] (a 3d frustrated Ising model 115411 for 
the (2+l)d case) or to a field-theoretical model similar to those studied in ref ]45l it is possible to 
show that for small K X y, valence-bond-solid states are realized in the background of polarized 
spins moments (m per site), as depicted in Fig sj3] and [4] (for (1+1) an d (2+1) dimensions respec- 
tively). These states support kink-like excitations with the fractional quantum number 8S Z = 1/2. 
The Id version of this valence-bond-solid plateau has been observed numerically ||55|]. It is inter- 
esting to note the apparent similarity of the Z 4 -vortex structure shown in Fig@]with those studied 
recently by other authors 1 56. Svfl in the context of fractionalized excitations in (2+l)d electron/spin 
models. In addition to these exotic states, more conventional 'collinear' (i.e. CDW) phases may 
appear when off-diagonal parts of M?? are sufficiently large. 

Since the 3d frustrated Ising model (equivalent to the (2+l)d version of eq. (|45l) with K X y = 0) 
has an order-disorder transition [|54|l . the phase diagram in the (2+l)d case is richer than its (l+l)d 
counterpart; a sufficiently small g m places the frustrated Ising model in a high-temperature param- 
agnetic phase (in this language the ordered phase with broken square-lattice symmetry translates 
into the columnar valence-bond- solid phase) which may be identified with the fractionalized phase 
of the model (I44al44bl) with the full space-group symmetry. This suggests that non-crystalline 
plateau states are possible in d > 2 while in Id, as implied by the LSM arguments, the appearance 
of plateaus is closely tied to the formation of crystalline states. (The dual vortex picture affords 
us with an alternative, simple description of the fractionalized excitations (spinons) which char- 
acterize this phase: they are the 2n vortices of the dual Ginzburg-Landau-type field theory which 
describes the condensate of doubly quantized vortices[45]. The Z 4 -kinks in the valence-bond- soid 
phase, while requiring a subtler analysis, can also be treated within a similar field theory language. 
) We summarize the above discussion in the form of a phase diagram, depicted schematically in 
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FigJH 

In concluding this section, we note that the two cases considered above, S — m E Z and 
S — m E Z + 1/2, each fall into the class of effective theories for easy plane antiferromagnets 
(with no magnetic fields and hence m — 0) with integer and half-integer !^ S. As we have 
already mentioned, however the present problem allows us to further extend this approach to other 
values of S — m. For instance the argument applies after appropriate modifications to the case with 
S — m = Z + 1/3; now the original boson is coupled to a Z 3 gauge field, and is fractionalized in 
such a way that it carries the 'charge' 5S Z = 1/3. The detailed analysis for general rational values 
of S — m will be reported elsewhere. 



^IGT 








FIG. 5: Schematic phase diagram of the fractionalized boson coupled to Z2 -gauge theory for (1+1 )d (inset) 
and (2+l)d. The large-KxY region is dominated by an XY-ordered (or, superfluid(SF)) phase, which, in 
(l+l)d, is replaced with an algebraic Tomonaga-Luttinger(TL) liquid phase. In the small-i^xY plateaus 
region, there is always a valence-bond(VB) phase with broken translation symmetry. Both in (l+l)d and in 
(2+l)d, the VB phase supports topological excitations (kinks in (l+l)d and Z4 vortices in (2+l)d). On top 
of this, there is & featureless plateau phase with fractionalized (i.e. 5S Z = 1/2) bosons in (2+l)d. 



VI. CONCLUSIONS 

The message of this article is twofold: (1) geometric phases provide a natural language to un- 
derstand how commensurability conditions impose constraints on the generation of magnetization 
in antiferromagnets, and (2) the analysis based on this language suggests a phase structure which 
accommodates a wealth of exotic phases. 

In the earlier portions of the article we found that a geometrical interpretation of the LSM 
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approach emerges when we view the spatial modulation of the slow variables as an " adiabatic 
evolution" along the spin chain. Continuing to focus on geometric phases, we were subsequently 
lead to effective field theories featuring vortex Berry phase factors which are valid in arbitrary 
dimensions. Taken together, these findings afford us with an explanation of a curious aspect of the 
LSM-style momentum-counting argumentfl which was mentioned back in subsection IIII At 
its ability to produce the correct criteria for the existence/nonexistence of a magnetization plateau 
despite the seemingly remote connection between the plateaus and the momenta of charge-neutral 
(SS Z = 0) excitations. The resolution lies in the fact that the LSM momentum shift and the 
action 5gp rtex governing the low-energy behavior via topological interference stems from a single 
quantity-the Berry phase of the canted spins. 

The approach exploited in the preceeding two sections is intended to serve as the basis for 
seeking as-yet undetected novel phases whose existence our theory suggests. Indeed the Z 2 gauge 
theory of section |Vl which is closely linked to the issue of fractionalization en route the pairing 
44J] strongly suggests that antiferromagnets in an external magnetic field provides 



of vortices fl43 



a promising playground to search for fractionalized spin liquids in spatial dimensions higher than 
one. (A magnetization plateau at m = 1/2 in translationally invariant S = 1 Hamiltonians on 
a 2d square lattice, a signature of such a phase, may in this regards be an interesting problem to 
pursue numerically as well as experimentally.) How the perturbation brought on by spin currents 
as briefly discussed towards the end of section HVl affects the magnetization property is another 
problem that we believe warrant further investigations. We hope to return to such issues in the 
future. 
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APPENDIX A: MORE ON THE BERRY PHASE APPROACH TO THE OYA ARGUMENT 

In this appendix we discuss how the arguments of section [III] are to be extended when one 
wishes to deal with a superposition of spin coherent states |{n(s)}). While the case treated in 
the main text already captures the essential geometric properties inherent in the LSM approach, 
this generalization intends to fill in the gap between the treatment of the actual ground state of 
a Heisenberg antiferromagnet (which possesses rotational invariance) with that given in the main 
text for a spin coherent state which is, in a semiclassical sense polarized in a particular orientation. 

Let \^f) then be our ground state. We can generally expand this state in the (overcomplete) spin 
coherent states basis where in accordance with the main text, the low energy physics is 

described in terms of the angular field 4>(x) specifying the in-plane orientation of the staggered 
spin component perpendicular to the magnetic field. 



The wave functional = ({(j)(x)}\^) is subject to the normalization condition 

/ T>(j)(x)\^[(j)(x)}\ 2 = I. We will adopt the action <S e// [0(r, x)] given in the text (eq.([23])), though 
the arguments to follow apply to a wider variety of systems (as long as the coefficient of the 
topological term remains the same). The Hamiltonian derived from this action is 



where k{x) = —ijj^j is canonically conjugate to 4>(x). The appearance of the gauge field-like 
piece S — m in the kinetic energy term has its origin in the topological term. In the superfluid 
analogy, it is simply the offset value of the superfluid density. The functional Schrodinger equation 
reads H{tt{x), <f>(x))V[<j>(x)] = E^l[<j>{x)]. 

As before we wish to compare the expectation values of the crystal momentum in the ground 
state and the LSM-twisted state. We recall that for the pure spin coherent state treated in sec- 
tionUm it was for this purpose essential that the quantity ({n(x)}|T|{n(x)}) pick a Berry phase 




(Al) 




(A2) 
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factor associated with a round trip experienced by the vector n. It turns out however that here 
the corresponding expression, (^\T\^) = J Vcj)(x)^*[(f)(x — a)] \l/ [</>(#)] is not equipped with 
the anholonomy one expects to find when <p is sent on a round excursion. This can be traced to 
the single- valuedness of \!/[</>(a;)] which implies that the wave functional be invariant under a 2%- 
shift of <j){x). Hence as it stands one cannot straightforwardly incorporate the Berry phase effect, 
which we know from the main text to be crucial in obtaining the correct expression for the crystal 
momentum. 

The resolution comes from seeking an analogy to the quantum mechanical problem of a particle 
on a ring threaded by an Aharonov-Bohm (AB) flux[|58ll to which our problem reduces when K± = 
0. As in the treatment of the latter, the AB-like anholonomy can encoded into the wavefunctional 
through a (singular) gauge transformation 

\j> [0( x )] = e ~* I dx{s-mWx) ^ (A3) 



which eliminates the AB-gauge field from the Hamiltonian, i.e. 



n 



dx 



1 .2 , ^l/o i\2 



(A4) 



at the price of introducing a twisted boundary condition for the transformed wavefunctional (such 
that a phase change of — 2n(S — m) is experienced as <p(x) changes by 2tv). 
We thus deduce the appropriate generalization of ({n(x)}|T|{n(x)}) to be 



(tf|T|tf) = / V(t)(x)^* 



a)]%(a;)] 



(A5) 



V(f){x)e-^ dx{s - m)d ^*[ct){x - a)]tf [0(z)], 

where we resorted to continuum notations in the final expression. Notice that the phase factor 
in the second line reproduces the correct Berry phase while the part (^f[(j)(x + a)])*\I/[0(x)] is 
single- valued. 

Introducing the unitary operator £/lsm = exp[— f dxx7r(x)} and using the relation 
uI sm TUlsm = exp[— J dx7t(x)]T, it is easy to see that the counterpart of eq. (|A5l ) for the 



LSM-twisted state reads 



(*I^sm^lsm|*> 
V(j)(x)& 

V(j)(x)ty* 



(x — a) 



2tt 

T 



( ^ 271 
[x — a) —a 



*[0(x)] 



(A6) 



^ e i2n(S—m)^—i J dx(S—m)d x (f> 
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Due to its single- valued nature, \!/[</>(x)] basically factorizes into factors of the form 
Yjn(x)& ( ^n{x)e in ^ x ^ x ' > coming from each point x. The effect of shifting the angular field cfi(x) 
by ^-a at each of the L/a sites thus merely amounts to a net factor of unity, i.e. ty[<f>(x) + ^a] = 
Hence we have 

(^I^sm^lsmI^) 

= J V<f)(x)y*[(f)(x - a))^[<f)(x)] (A7) 

x ^i2ix(S-m) —i f dx(S-m)d x (j> 

A comparison of eqs. (IA5l) and (IA7I) yields the desired quantization rule. 



APPENDIX B: BOSON- VORTEX DUALITY IN THE PRESENCE OF A TOPOLOGICAL 
TERM 

Here we provide the main steps leadings to the Berry phases of eq. (!26cl) . which makes extensive 
use of the boson-vortex duality transformation (performed at the first quantization level). Though 
the technique is standard, we feel that it may be worthwhile to illustrate how the presence of a 
topological term brings about modifications and gives rise to vortex Berry phases which ultimately 
dominate the physics at low energy. We also sketch how to retain this Berry phase effect when one 
switches from the first to the second quantized description of the vortices. 



1. First quantization 

For simplicity we work in the continuum formulation while being careful to take proper account 
of the Berry phases; carrying out a lattice duality transformation will of course produce the same 
results. We begin by introducing a Hubbard-Stratonovich auxiliary vector field J M = ( J T , J) and 
rewrite the imaginary-time lagrangian density of the effective theory <S[0(r, r)] = J drd d rC as 

C —i I J T H -j— d T 4> + 



+ iJ-V(f) + 



a d ) " ' T ' 2K T 
J 



(Bl) 



2K ± 

This is followed by a decomposition of the phase field into components with and without vor- 
ticity, 4> = V + 4> r , where (d^d u — d^d u )(fi Y ^ 0. The vorticity-free component <fi r can be safely 
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integrated out, which yields the constraint 

<V M = 0, (B2) 
where J /i = J fi + S^^r, and the lagrangian density becomes 



2K r 2K\ 



C = U^fa + v ~ or /' ' + —. (B3) 



J, 



The divergence-free condition, eq. (|B2|) . is solved explicitly in terms of a new field whose form 
depends on the dimensionality: 

e^ u d v <p (d = l) 

e^xd v b x (d = 2) (B4) 
^tMuXpd u bxp (d = 3) 

where ip, b^, and 6 M „(= —b u ^) are all vorticity-free. Here on the procedures will be described 
separately for each dimensionality d. 

(i) d=\: We define the density of spacetime vortices p v via (d T d x — d x d T )(f) v = J2j2nqj5(T — 
rJ)S(x — Xj) = 2np v , where (rj,xj) is the spacetime coordinate of the j-th vortex event (the 
notation Xj which appeared in eq. (|26cl) in the main text corresponds here to a; J. ). After an 
integration by parts and a shifting of the t^-field, ip = cp + ^j^x, we arrive at 

p + — - — xj. (B5) 

Integrating over <p we obtain 

C = 7r 2 p v —p v + t2TT P v ^^x. (B6) 
— o 2 a 

The first term on the right hand side is the intervortex Coulombic interaction for which the kernel 
l/(— d 2 ) was defined in the main text. The second term represents the vortex Berry phase terms, 

v S—m _v 



X-. 



i.e. 12tx J drdxp v ^x = ^vr^g^ 
(ii) d=2: We start by establishing some notations. The 3-vector b^ in eq. (|B4|) is often referred 
to as the vortex gauge field or the dual gauge field (note the invariance of eq. (IB4l) (in (2+l)d) 
with respect to the gauge transformation b^ — > b^ + d^x). Dual electric/magnetic fields can be 
constructed from these gauge fields: 

E dual = (jB dual ; ^dual) = ^ _ dxb ^ drby _ d 

(B7) 
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They are related to the auxiliary field through 

J T = B dua \ 

(B8) 

j pdual t pdual 

J x — i J y — 

Upon introducing the vortex 3-current = 7^re M „ A ct<9 A v which couples to the dual gauge 
field, the lagrangian density can be rewritten as 

B ^_S-m\ + 1 (Edual)2+ . 2 r (B9) 



2K T V a 2 J 2K ± y ' 11 M 

We then divide 6 M (choosing a suitable gauge) into two parts, viz. = b^+db^, where b^ generates 
a background magnetic field d x b y — d y b x = while Sb^ represents the deviation from this offset 
value, 5B dual = B dual — = d x 5b y — d y 5b x . The gauge is so chosen that b^ does not modify 
the dual electric field E dual (a simple choice being (b T , b x , b y ) = (0, 0, ^?rx)). Having set up the 
necessary notations, we can recast eq. (|B9l ) as 

C = ^^^^ + ^^^^ + l2n6b ^ + i2 ^ J l- (B10) 
On integrating over 5b^ we are left with 

C = /^coulomb + ilnb^r^ (Bl 1) 

where the first term on the right hand side is again the Coulombic interaction (now between seg- 
ments of the spacetime vortex loops) which is of the form n 2 ^(rp)^ — ^^)Jl- The second 
term assigns a Berry phase to each spacetime vortex loop, as we will now see. Since the vortex 
current field is nonvanishing only on such loops (with a magnitude equal to the vorticity qj), 
we can convert the spacetime integration of this term into a sum of contour integrals along each 
vortex loop (we use below a short-hand notation V to denote a (d+l)-vector V^) : 

^° P rtcx = i2w J drdhb^l 

The parameter s(e [0, 1]) is used to specify the distance along each loop contour Cj. Using Stokes' 
theorem this contour integral can in turn be expressed as an area integral. In vectorial notations 
i2n J2j Qj §c d>Xj'b(xj) = i2ir J2j l] f D dA-rotb where dDj = Cj and dA is the area element on 
the two-dimensional domain Dj. Recalling that by definition rot b = ^?re T , the final expression 
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for the Berry phase term is i2ir(S — m) ^ . q]A* y , where A xy is the area bounded by the xy-plane 
projection of the loop Cj. 

(iii) d = 3: The procedures involved are natural extensions of the d = 2 case. (Readers seeking 



more details on the framework specific to d=3 can consult the work of Zee[36].) 
The lagrangian density corresponding to eq. (|B9l) is 

\ / g _ m \ 2 i 
£ = {tijkdibjk ^— J + j-^r- [(e xa p 7 d a bp 7 ) 2 + (e^^fr^) 2 ] +i27rfe Atl/ J^, (B13) 

where the suffixes in roman letters stand for spatial indices, and information on the vortex 
world sheets is encoded in the 2-form T 1 = -^e^xpdxdp^. Within a worldsheet parametriza- 
tion scheme using a set of two parameters (s, u), the latter can also be written as J* = 
^jQ] J ds J due^daxj^db xJ^S^ (x — Xj(s,u)), where a, b axe s or u (compare with expres- 
sion for vortex current employed in eq. (IB12l) ). Since each worldsheet is nothing but the spacetime 
trajectory of a vortex loop living in 3d space, JX, can be regarded as a vortex loop current. As in 
the 2+ Id case we make the decomposition b^ = b^ u + 5b ^ with the background field satisfying 

777, 

eijkdib jk = — , (B14) 

(X 

together with the condition that e^ v xd^b u \ = when one of the indices is r. Integrating out 5b^ v , 
we get 

£ = ^Coulomb + i ZnJ^J)^. (B15) 

The first term is the current-current interaction between worldsheets, while the second is the Berry 
phase term, which contributes to the action a sum of surface integrals on each worldsheet. The 
latter converts, upon an application of the Stokes-Gauss theorem into 3-volume integrals of the 
flux e^\pdj)\ p . Since the only nonzero flux component comes from eq. (|B14l) . we find that the 
vortex Berry phase amounts to i2ir(S — m) <?J^7^' where V* yz is the 3-volume of the interior 
of the j-th worldsheet projected onto the rq/z-subspace. 

2. Second quantization (vortex field theory) 

The dual theories derived above can be used to obtain the vortex field theories of eq.(l27l). We 
give a shortcut version of this procedure for the case of d — 2 which, as mentioned in the main 



text, is interesting in view of its relation to the work of Lee and Shankar[42] (the steps for the 



cases d = 1 and d = 3 are essentially the same). We first make the following relabelling in 
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eq. dBlOD : 5b^ — > b^, 5B dual — ► B, and E dual — > E. The action is then decomposed into two parts, 
S = 5kin + Sj.b, where 5kin = / d 3 x[j^-B 2 + is the Maxwellian kinetic energy term, and 

Sj.b gives the coupling between the vortex current and the dual gauge field. We write the latter as 

Sj. b = -i2ii J d 3 xJ^(d li (p — b f j, — 6 l j,ypx), (B16) 

where as before the "superfluid density" in (2+l)d is p = ^z 1 , and we have introduced a new 
scalar field tp whose role is to impose the continuity of the vortex current[21], d^J^ = 0. We now 
place our system on a spacetime lattice. In so doing, note that the vortex current can be written 
(as in eq. (|B12l) ) as = ^ gj § c (x — #J(s)). In the discretized lattice version of 

the theory, each segment of the vortex current can be represented by an integer-valued 3-vector 
= <l]~ctr defined on the dual lattice. 

The partition function of a grand-canonical ensemble of spacetime vortex loops is thus 

Z= [{Jl dtp n db ntl ) II e-^^e - *^^^- 1 ^-'^^-^, (B17) 

n,n {t n ^} n,fi 

where n = (n T ,n x ,n y ) is the site index, and A the line tension of each vortex loop. Using the 
Poission resummation formula, we can trade the sum over t nfl with an integration over a continuous 
variable X nfJr , at the price of introducing another set of integers {m n ^}. After the X nil -integration 
is carried out, we obtain 

Z = J (Y[dfi n)1 db nfl ) J2 J^e-^^e-t^ 2 ^^- 5 ™^^"^] 2 , (B18) 

nfi {rn nfl } rifi 

which is just the Villain form of 

Z = I ([[dfi nf tdb n ^ JJ e -^[M e -^co S [2 7 r(A^-6 M -«5 MH pn,)]_ (B19) 
nfi nfi 
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